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Abstract 

We investigate by means of Monte Carlo simulation and Finite-Size Scaling anal- 
ysis the critical properties of the three dimensional 0(5) non linear a model and 
of the antiferromagnetic RP 2 model, both of them regularized on a lattice. High 
accuracy estimates are obtained for the critical exponents, universal dimensionless 
quantities and critical couplings. It is concluded that both models belong to the 
same Universality Class, provide that rather non standard identifications are made 
for the momentum-space propagator of the RP 2 model. We have also investigated 
the phase diagram of the RP 2 model extended by a second-neighbor interaction. A 
rich phase diagram is found, where most phase transitions are first order. 
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1 Introduction 



Universality is sometimes expressed in a somehow defectively simple way: some 
critical properties (the universal ones) of a system are given by space dimen- 
sionality and the local properties (i.e. near the identity element) of the coset 
space Q /TC, where Q is the symmetry group of the high temperature phase 
and TC is the remaining symmetry group of the broken phase (low tempera- 
ture). As we shall discuss, the subtle point making the above statement not 
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straightforward to use, is that Q needs not to be the symmetry group of the 
microscopic Hamiltonian, but that of the coarse-grained fixed-point action. 

On the spirit of the above statement, some time ago [1] a seemingly complete 
classification was obtained of the Universality Classes of three dimensional 
systems where £=0(3). In this picture, a phase transition of a vector model, 
with 0(3) global symmetry and with an 0(2) low temperature phase symme- 
try, in three dimensions must belong to the 0(3)/0(2) scheme of symmetry 
breaking (classical Heisenberg model). In addition, if 7i=0(l)=Z 2 is the re- 
maining symmetry, the corresponding scheme should be 0(4)/0(3) which is 
locally isomorphic to 0(3)/0(l). 1 This classification has been challenged by 
the chiral models[2]. However, the situation is still hotly debated: some au- 
thors believe that the chiral transitions are weakly first-order [3], while others 
claim[4] that the Chiral Universality Class exists, implying the relevance of 
the global properties of Q /H. 

In this letter, we shall consider the three dimensional antiferromagnetic (AFM) 
RP 2 model [5,6,7,8], a model displaying a second-order phase transition and 
escaping from the previously expressed paradigm. It is worth recalling [9,10] 
that one of the phase transitions found in models for colossal magnetorresis- 
tance oxides [11] belongs to the Universality class of the AFM RP 2 model. 
The microscopic Hamiltonian of this model has a global 0(3) symmetry group, 
while the low temperature phase has, at least, a remaining 0(2) symmetry [9]. 
We will show here that the model belongs to the Universality Class of the 
three dimensional 0(5) Non Linear a model. Some ground for this arises from 
a hand- waving argument, suggested to us by one of the referees of Ref. [9] (see 
below) . 

The Universality Class of the three dimensional 0(5) Non Linear a model 
has received less attention that 0(N) models with N = 0,1,2,3 and 4. In 
spite of that, it has been recently argued that 0(5) could be relevant for 
the High-temperature superconducting cuprates [12]. Nevertheless, pertur- 
bative field theoretic-methods have been used to estimate the critical expo- 
nents [13,14,15,16]. From the numerical side, only a rather unconvincing Monte 
Carlo simulation [17] was available until very recently. Fortunately, there has 
been a recent, much more careful study [18]. Yet, the scope of Ref. [18] was to 
determine whether an interaction explicitly degrading the 0(5) symmetry to 
an 0(3) © 0(2) group was relevant in the Renormalization-Group sense. To 
that end, those authors concentrated in producing extremely accurate data on 
small lattices. 

Our purpose is to study in greater detail the critical properties of the three di- 
mensional 0(5) Non Linear a model, and of the AFM RP 2 model. We improve 

1 This statement assumes that the global properties of the coset G/TC are irrelevant, 
only the local properties matter. 
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over previous studies of both models, obtaining more accurate estimates for 
critical exponents, universal dimensionless quantities and non universal criti- 
cal couplings. As symmetries play such a prominent role, we will also explore 
the possibilities of changing those of the low-temperature phase by adding a 
second-neighbors coupling to the Hamiltonian of the AFM RP 2 model. 



2 The models 



We are considering a system of A^-component normalized spins {vl} placed 
in a three dimensional simple cubic lattice of size L with periodic boundary 
conditions. The actions of our lattice systems are 

<i,j> <i,j> 



where the sums are extended to all pairs of nearest neighbors. Our sign con- 
vention is fixed by the partition function: 

Z=[l[dv l e- S , (2) 



dv being the rotationally invariant measure over the A^-dimensional unit sphere. 

To construct observables, in addition to the vector field vl, we consider the 
(traceless) tensorial field 

rf = v?v?-^5 a ^ a, P — 1, . . . , N . (3) 



The interesting quantities related with the order parameters can be con- 
structed in terms of the Fourier transforms of the fields (fi = Vi, r^) 

Kp) = ^Y.^h- (4) 



For RP^ -1 models, the local gauge invariance Vi — > — implies that the 
relevant observables are constructed in terms of the tensor field. However, for 
O^) we have found very interesting to consider as well quantities related 
with the tensor field. 
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We construct the scalars (under global O(N) transformations and spatial 
translations) 

S y (p) = %{p) ■ %\p), S T (p) = tr r(p)r*(p) , (5) 



which, in addition to the action, are the only quantities measured during the 
simulation. Their mean values yield the propagators: 

G T ,v(p) = L*(Sr,v(p)) ■ (6) 



In the thermodynamic limit and at the critical point, the propagator is ex- 
pected to have poles at p { = (0, 0, 0) and, for the antiferromagnetic model, at 

p s = (7T,7T,7r): 2 



where £||5p|| ^ 1> an d the exponents rf and rf correspond to independent 
wave function renormalization at each pole. Note that close to the critical 
point £ f and £ s are expected to remain proportional to each other (this will 
be explicitly checked numerically). 

The (non-connected) susceptibilities are simply: 

X = G(p ) • (8) 



In a finite lattice an extremely useful definition of the correlation length can 
be obtained from the (discrete) derivative of G(p). Using 5p = (2ir/L, 0,0) 
one obtains[19,20] 

_ ( G(p )/G( Po + 5p)-l \ 1/2 
^"v. 4sin 2 (7r/L) ) ■ W 



We also compute the cumulants 

u. = ^ ■ (10) 



2 In the remaining part of this section, if a subindex V or T does not explicitly 
appears, it will imply that the equation is valid both for vector or tensor quan- 
tities. The same convention will apply for superscript (f) (ferromagnetic) and (s) 
(staggered). Staggered quantities are useful only for antiferromagnetic models. 



4 



Finally, the energy per link is 



e = (s/(-m 3 )) . 



(ii) 



We have computed /^-derivatives of observables through their connected cor- 
relation with the action. Furthermore, we have extrapolated mean-values from 
the simulation coupling, to a neighboring value of f3 using the standard reweight- 
ing techniques, that cover all the relevant part of the critical region [20] . 

The relationship between the 0(5) model and the RP 2 model arises from the 
Landau- Wilson-Fisher Hamiltonian for the RP 2 system [9]. Indeed, at the 
Mean Field level [9,6], the ferromagnetic quantities are simple functions of the 
staggered ones. This suggests to construct the Landau- Wilson-Fisher Hamil- 
tonian from the staggered magnetization, which is a traceless, real, symmetric 
3x3 matrix: 



Note that M s has 5 independent quantities. It is therefore a simple matter to 
obtain a five-components real vector v such that v 2 = trAi 2 . The less trivial 
part regards the fourth-order interaction terms. In principle, the 0(3) symme- 
try of the microscopic Hamiltonian would allow for a trAij term and a [tr.M 2 ] 2 . 
Surprisingly enough, both terms are proportional to iy 2 ) 2 . Thus, assuming 
that sixth-order terms are irrelevant, the Landau- Wilson-Fisher Hamiltonian 
is expected to have a 0(5) symmetry group and both models belong to the 
same Universality Class. This does not only implies that both models have 
the same critical exponents but also that the L — > oo limit of U^^p and of 

U% RF (evaluated at their respective critical couplings) coincide. 



3 Numerical methods 

In the 0(5) model we have studied lattice sizes L — 6, 8, 12, 16, 24, 32, 48, 64, 96 
and L = 128, at f3 = 1.1812. We have combined a Wolff's single cluster 
update with Metropolis. Our elementary Monte Carlo step (EMCS) consists 
of (10L + 1) Wolff's cluster update and then a full-lattice Metropolis sweep. We 
take measurements after every EMCS. Since the average size of clusters grows 
as L 2 ~ n pa L 2 , 80% of simulation time we are tracing clusters for all L, while 
Metropolis accounts for 10% of time and measurements for the remaining 10%. 
The total simulation time time has been the equivalent to 600 days of Pentium 
IV at 3.2 GHz. The number of EMCS ranges from 10 8 for L = 6 to 1.4 x 10 6 



M?" = £ (-i) 



Xi+yi+Zi a,/3 



(12) 
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at L — 128. The integrated autocorrelation times for the susceptibility and 
for the energy are smaller than 1 EMCS for all the simulated lattice sizes. 

Since we are interested in high accuracy estimates, we have used double pre- 
cision arithmetics. One also needs to worry about the pseudo random-number 
generator. We have therefore implemented a Schwinger-Dyson test. It turned 
out that the 32-bits Parisi-Rapuano pseudo random- number generator [21] 
produces biased results. Either the Parisi-Rapuano plus congruential genera- 
tor [22] or the 64 bits Parisi-Rapuano generator cured this bias. The 64 bits 
Parisi-Rapuano generator is faster and it has been our final choice. 

For the antiferromagnetic RP 2 model, no efficient cluster method is available. 
We have simulated in lattice sizes from L = 8,12,16,24,32,48 and L = 64 
at P — —2.41. We used a multi-hit Metropolis sequential algorithm. Making 
a new spin proposal completely independent from the previous spin value, we 
achieve an acceptance of about 30%. We have used 2 hits what ensures a 50% 
acceptance. The observables have been measured every two Metropolis full 
lattice-sweep (our EMCS). 

The number of EMCS ranges from 10 8 for L = 8 to 7 x 10 8 for L = 64. In 
units of the integrated autocorrelation time r (for the order parameter) we 
have more than 10 6 r for L = 64. The data up to L = 48 were obtained in 
Pentium IV clusters (simulation time was roughly equivalent to a 1000 days 
of a single processor). For the largest lattice, data were obtained in the Mare 
Nostrum computer of the Barcelona Supercomputing Center (simulation time 
was roughly equivalent to 3000 days of a single processor). 

We perform a Finite-Size Scaling analysis, using the quotients method [5,6,20]. 
In this approach, one compares the mean value of an observable, O, in two 
systems of sizes L\ and at the value of (3 where the correlation- length in 
units of the lattice sizes coincides for both systems. If, for the infinite volume 
system, {0)((3) oc \(3 — f3 c \~ x ° , the basic equation of the quotient method is 

= (-r) (i+AoLr + ...), (13) 

«(il,/3) H 

where the dots stand for higher order scaling corrections, v is the correlation- 
length critical exponent, — u is the (universal) first irrelevant critical exponent, 
while Ao is a non universal amplitude. In a typical application, one fixes 
the ratio s = L 2 /Li to 2, and consider pairs of lattices L and 2L. A linear 
extrapolation in L~ w is used to extract the infinite volume limit. One just needs 
to make sure that the minimum lattice size included in the extrapolation is 
large enough to safely neglect the higher-order corrections. Of course, any 
quantity scaling like £ at the critical point, such as LU4, may play the same 
role in Eq.(13). However, usually £ yields smaller scaling corrections than U4. 



MM = (°(^ L 2)) 

~ WM)) 
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The extrapolation method based on Eq.(13) is feasible for the antiferromag- 
netic RP 2 model. Unfortunately, for the 0(5) model the amplitude Aq is sur- 
prisingly small. In fact, resummation of the e-expansion yields oj = 0.79(2) [13], 
while blind use of Eq.(13) on our numerical data would predict u « 2. We 
have then considered an additional correction term, AoL~ a . The exponent 
a is an effective way of taking into account a variety of higher-order scaling 
corrections of similar magnitude (an L~ 2uJ contribution, subleading universal 
irrelevant critical corrections, analytic corrections, effects of the non-linearity 
of the scaling fields, etc. [20]). Its utility will be in that it allow us to give sen- 
sible error estimates for the infinite-volume extrapolations, instead of bluntly 
taking A = 0. 

The most precise way of extracting the critical exponent, u, and the critical 
point j3 c is to consider the crossing point of dimensionless quantities such as 
£/L and U±. Indeed, comparing their values in lattices L\ and L 2 , one finds 
that they take a common value at 

' Jc ~' Jc + B (L 2 /L 1 y/^-i 1 + '--' (14) 

The non universal amplitude B depends on the considered dimensionless quan- 
tity. Again, one usually take pairs of lattices L and 2L, and extrapolates to 
infinite volume using Eq.(14), maybe performing a joint fit for the crossing 
points of several dimensionless quantities. Again, for the 0(5) model the am- 
plitudes B are exceedingly small, and we need to add to Eq.(14) an analogous 
higher-order term, where a plays the role of u, and with amplitude B. 



4 Results for the 0(5) model 



The first step is the location of the critical point and the scaling corrections 
exponent. In table 1 we show the crossing points f3^ ,2L for the dimensionless 
quantities £y/L and U^y. To study the finite size corrections, we need to fit 
them to 

Pc' 2L « & + BL~^ l/u + BL- a - l/u . (15) 



Fixing B = yields uo larger than 2 which is unacceptable given the field 
theory estimate uj = 0.79(2) [13]. Our interpretation is that B is too small to be 
observed even with our 6-digit accuracy. We have therefore fixed 00 = 0.79(4) 
and we have taken as fit parameters /3 C , B, B and a. We have doubled the 
field theory error in u for safety. To further constrain f3 c (and a), we have 
performed a joint fit of the crossing points for both £y/L and U^y, with the 
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1.181365(3) 


0.28146(3) 


1.181371(4) 


1.06976(2) 


48 


1.181366(3) 


0.28155(5) 


1.181373(5) 


1.06972(3) 


64 


1.181362(4) 


0.28142(9) 


1.181370(6) 


1.06982(6) 



Table 1 

Effective critical points, /3c ' 2L , of the 0(5) model, from the crossing points of the 
dimensionless quantities £v/L and U^y obtained in lattices of sizes L and 2L. We 
also show the value of both quantities at their respective crossing point. 

same f3 c . For this model we always fit for L > 8 . The results are (in all fits 
reported in this work, the full covariance matrix was used): 

P c = 1.1813654(19) , a = 2.21(17) , x7 d o f = 7.1/8 . (16) 



Notice that, for using equation (15), an estimate of v is needed. Fortunately 
a rough estimate v ~ 0.78 (see below) is enough, given the uncertainty in oo 
and a. As for the amplitudes of the leading correction term, we find 

5 5v/L = 0.004(6) , B U4V = -0.001(3) , (17) 



while the amplitudes B are of order one. We then see that the L a term is 
crucial in order to obtain a sensible error estimate in the L — > oo extrapolation. 

At this point we may obtain two universal quantities, £7| v and £y/L, namely 
the L — > oo limit of U±y and £y/L evaluated exactly at the critical coupling. 
Again, due to the smallness of the leading scaling-corrections, we extrapolated 
to L — > oo using the following functional forms: 

U 4y (^' 2L , L) » Uly + C Uiy L-" + C UtiV L-° , (18) 



and 



j; + Qv/l j 



+ c t 



(19) 
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Our numerical estimates for U^\(f3% ,2i , L) and ^y((3^' 2L , L)/L are displayed 
in the third and fourth columns of table 1. Although scaling-corrections are 
tiny, they can be clearly observed. Using Eqs. (18,19) we obtain 



| = 0.28145(13) = ^), W AN = 1.06978(5) (^ = M) , (20) 
while the amplitudes of the leading scaling-corrections are 

Q v/L = 0.002(3) , C U4tV = 0.0002(6) . (21) 



L 




77v 


77T 


6 


0.7963(3) 


0.04343(6) 


1.3499(3) 


8 


0.7894(4) 


0.03837(7) 


1.34254(12) 


12 


0.7849(4) 


0.03554(5) 


1.33780(10) 


16 


0.7833(4) 


0.03462(6) 


1.33589(11) 


24 


0.7819(7) 


0.03430(10) 


1.33455(19) 


32 


0.7802(14) 


0.03396(16) 


1.3332(2) 


48 


0.7817(19) 


0.0339(2) 


1.3322(4) 


64 


0.781(4) 


0.0341(4) 


1.3321(7) 



Table 2 

L dependent effective values of exponents u,r}y and t\t for the 0(5) model, calculated 
from Qq 2L . 

To obtain the critical exponents, we consider the operators <9g£v and xv,t 
whose associated exponents are xg^ v — v + 1 and x Xv T = 7v,t = ^(2 — 
t^v.t) • Taking the base 2 logarithm of the quotients, see Eq.(13), we obtain 
the effective size dependent exponents shown in table 2. In order to obtain 
their infinite volume value, we use (13), including an explicit L~ a term in the 
fit: 

qL,2L = ro /u + AqL -u + ~ AqL -* ^ 



(note that we have absorbed a constant factor 2 X °/ U into the amplitudes A for 
scaling-corrections). We obtain 

v = 0.780(2), a = 2.15(19), A 9 ^ v = -0.04(7), x 2 /dof = 8.4/8 

7] Y = 0.03405(3), a = 2.27(19), A Xv = 0.0012(14), x7 dof = 8.5/8 (23) 

77t = 1.3307(5), a = 2.24(19), A XT = -0.0053(12), x 2 / dof = 9 - 6 / 8 
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Fig. 1. Size-dependent estimators for the anomalous dimensions of the 0(5) model 
as obtained from (L,2L) pairs, versus . We plot estimates from the crossing 
points of U±y an d of £v/L. We use the value u = 0.79, from Ref.[13]. 

Less than a 5% of the total error is due to the error in u = 0.79(4) for both v 
and tjy. For tjt it is about a 30%. There are two points to be made about the 
extrapolation: 

• The 0(5) model is not an improved action [23] (in the sense of exactly- 
vanishing leading scaling-corrections), since A XT is clearly non-zero (this is 
also illustrated in Fig. 1). Had we not considered the tensorial operators, 
this would have been completely missed. 

• It is somehow disappointing to compare the accuracy of the effective expo- 
nent v in table 2 with the error in the extrapolation (23). Indeed, could we 
safely set Aq^ = 0, the final result would have been v = 0.7813(4). 



5 Results for the antiferromagnetic RP 2 model 

As we said before, qualitative arguments suggest that the antiferromagnetic 
RP 2 model belongs to the 0(5) Universality Class. Our aim is to make the 
most astringent possible test of this hypothesis, thus we perform here an up- 
date of a previous study [6] of the RP 2 critical quantities. We report here 
largely improved estimates for critical coupling and exponents. Furthermore, 
we give estimates for the dimensionless quantity C/|'*, that can be directly 
compared with the U^ v obtained for the 0(5) model. 
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Fig. 2. Crossing points of the different dimensionless quantities, for the (L, 2L) pairs, 
as a function of L~ w ~ x l v ' . Lines are fits to Eq.(14), constrained to yield common 
values of j3 c and oj. The oj used in the X axis correspond to the optimal value. Note 
that the minimum pair plotted is (12,24). 

In this case, the extrapolation to the infinite volume limit is more standard 
(see for instance [24]) than for the 0(5) model, because the amplitude of 
the leading scaling-corrections are much larger in most cases. To estimate 
oj and (3 C , we consider pair of lattices L and 2L, performing a joint fit to 
Eq.(14) of the crossing points of all four dimensionless quantities, imposing a 
common value of (3 C and oj (see Fig. 2). To control for systematic errors due to 
higher-order corrections, we follow the following procedure. We perform the 
fit using data for L > L min , seeking a value of L min where a reasonable value of 
X 2 /dof is found. Furthermore, we require that the fit performed for L > L min 
yield compatible results. In that case, we report the central value from the 
L > L m[n fit, but taking the enlarged errors from the L > L m[n fit. We found 
that L min = 12 is enough for the extrapolation of f3 c . We obtain: 

(3 C = -2.40899(13) , oj = 0.78(4) , xVdof = 8.5/10 . (24) 

Once we have determined oj, we proceed to extrapolate f/J*, and the critical 
exponents, using the analog of Eqs. (18, 19,22) without the effective L~ a term. 
Although one could consider all four types of crossing points, /5 ( f' 2L , the re- 
sulting quotients would be highly correlated, making join fits scarcely useful. 
We concentrate on the crossing point of £ S /L, which seems the most natural 
quantity, as we are dealing with an antiferromagnetic model and it can be 
obtained using only the two points correlation function. We have checked that 



other choices for (3^ ,2L yield compatible results, with slightly larger errors. Our 
extrapolations are shown together with the effective L dependent estimates in 
table 3. Error estimates in the extrapolation include the effect of the uncer- 
tainty in uj. For exponent v, scaling corrections are completely buried in the 
statistical errors. We extrapolated with a simple linear fit, using L min = 8. The 
situation is rather different for 77 s . For that exponent, enlarging L m j n system- 
atically increases the asymptotic estimate. On the other hand, a fit quadratic 
in L~ u yields a linear term compatible with zero. The linear extrapolation 
with L min = 16 is identical to the quadratic extrapolation from L min = 8. 
This is the result indicated in table 3. As for rf, we have rather strong leading 
scaling corrections. Indeed, a fit linear in yields basically identical results 
for L m m = 8 and L m i n = 12 (this is the result reported in table 3). Further- 
more, a fit quadratic in L~ w including all points, yielded rf = 1.331(5). The 
extrapolation for C/|'* is equally simple. 



L 


R L,2L 
Pc,i'/L 


W 


V 


rf 




8 


-2.39892(18) 


1.06810(8) 


0.7848(12) 


0.0390(3) 


1.4155(6) 


12 


-2.40487(13) 


1.06805(11) 


0.7871(16) 


0.0364(4) 


1.3902(7) 


16 


-2.40670(12) 


1.06849(14) 


0.782(2) 


0.0357(6) 


1.3776(11) 


24 


-2.40792(9) 


1.06883(19) 


0.779(3) 


0.0351(7) 


1.3655(13) 


32 


-2.40846(7) 


1.06868(18) 


0.783(3) 


0.0337(7) 


1.3572(14) 


00 


-2.40899(13) 


1.0691(5) 


0.780(4) 


0.032(2) 


1.328(4) 



Table 3 

Size dependent estimators for the critical coupling and several Universal quantities, 
as obtained from (L, 2L) pairs, in the RP 2 model. The last row correspond to the 
infinite volume extrapolations. 

The extrapolation for other scale-invariant quantities, without an obvious cor- 
respondent in the 0(5) model, are: 

i_ = 0.5379(17) , ^— = 0.2236(15) , U f 4 '* = 1.3114(6) . (25) 

Lj Lj 



6 Next nearest neighbors coupling 

A rather subtle question regards the symmetry of the low-temperature RP 2 
antiferromagnetic phase [6,9]. A way of investigating this problem is to study 
the enlarged action 

<i,3> <«J> 



12 



1 

0.5 


-0.5 
-1 

(N 
CO. 

-1.5 

-2 
-2.5 

-3 
-3.5 

-6 -5 -4 -3 -2 -1 

Pi 

Fig. 3. Phase diagram of the extended antiferromagnetic RP 2 model, Eq.(26). 
where an additional second-neighbors coupling is considered. 

The phase diagram for f3i < (Fig. 3) contains the following regions (spins 
are classified as even or odd, according to the parity of Xi + + z^): 

• PM: the usual (paramagnetic) disordered state, where the 0(3) symmetry 
of the action (26) is preserved. 

• 0(2): (say) even spins fluctuate almost parallel to (say) the Z axis, with 
random sense (local Zi symmetry), while odd spins fluctuate in the perpen- 
dicular plane (global 0(2) symmetry) 

• 0(1): two sublattices with ferromagnetic ordering in perpendicular direc- 
tions, with random sense (the local Z 2 symmetry, v { — > — vl is always pre- 
served) . 

• Skyrmion/Flux: the spins are parallel to the diagonals of the unit cube, 
so that they point out from/to the center (i.e. the propagator show three 
peaks at p = (tt, it, 0) and permutations). It is interesting to note the vec- 
torial version of this phase appear in models for colossal magnetorresistance 
oxides [25]. 

The most relevant results can be summarized as follows: 

• We have obtained critical exponents for several points along the PM-0(2) 
critical lines, with significantly less accuracy than for the 02 = model. No 
variation was observed within errors. 

• The 0(2)-0(l) critical line is repelled from the $2 = axis by the second- 
neighbors coupling. We interprete this as a competition with the order-from- 
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disorder mechanism [6,9] behind the PM-0(2) transition. 

• A naive analysis suggests that the 0(2)-0(l) transition line should belong 
to the XY universality class (consider first the limit f3\ = — oo, then the 
identity cos 2 9 = (1 + cos26>)/2 for the less ordered face-centered cubic 
sublattice). However, we have found that this transition line is first order, as 
revealed by the double-peaked histogram of the second neighbors energy. We 
are able to estimate a non-zero latent-heat up to f3\ = —6. At f3\ = —4 the 
double-peak structure is still easy to observe on small lattices. We presume 
that the whole line is first-order, although it could become very weak. 

• The skyrmion-PM transition lines turned out to be first-order at all the 
checked points. 

• Note that at f3\ = 0, we have two decoupled ferromagnetic RP 2 models on 
the face-centered cubic lattice (a model showing first-order transition, well 
known in the liquids-crystal context [26]). We should remark that a precise 
location of the triple point 0(2)-0(l)-PM is very difficult to achieve. 



7 Conclusions 

We have obtained high accuracy estimates of critical exponents and other 
universal quantities for the three-dimensional 0(5) and the antiferromagnetic 
RP 2 models, by means of Monte Carlo simulation, Finite-Size Scaling analysis 
and careful infinite-volume extrapolation. 

In the case of the 0(5) model the coupling to the leading irrelevant operator is 
rather weak, but non-vanishing, and one needs to consider higher-order scaling 
corrections to obtain sensible error estimates. In spite of that, our estimates for 
the critical coupling and the anomalous dimensions for the vector and tensor 
representations improve significantly over previous work (see table 4). 

For the RP 2 model, the leading scaling corrections are sizeable. It is amus- 
ing that we are able to obtain numerically (for the first time, we believe) an 
estimate for the (universal) scaling-correction exponent u = 0.78(4), of accu- 
racy comparable to the perturbative field-theoretical estimate[13] u = 0.79(2). 
As table 4 shows, within the achieved accuracy, both models seem to belong 
to the same Universality Class (for comparison, we also show results from 
the 0(4) Universality Class). To conclude this, one needs to accept that in 
the RP 2 model, the wavefunction renormalization for the propagator pole at 
p = (ir, 7r, 7r) is as for the the 0(5) fundamental field, while at p = (0, 0, 0) 
is as for the 0(5) tensor field. 

We have also obtained the phase-diagram of the RP 2 model extended with a 
second nearest-neighbors interaction. We have found a rich phase diagram. 
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Model 


Pc 


V 


V 


rf 


Ui 


RP 2 (this work) 


-2.40899(13) 


0.780(4) 


0.032(2) 


1.328(4) 


1.0691(5) 


0(5) (this work) 


+1.1813654(19) 


0.780(2) 


0.03405(3) 


1.3307(5) 


1.06978(5) 


0(5) (Rcf. [18]) 


+ 1.18138(3) 


0.779(3) 


0.034(1) 




1.069(1) 


0(5) (FT [13]) 




0.762(7) 


0.034(4) 






0(4) 


+0.935858(8)[27] 


0.749(2) [28] 


0.0365(10) [28] 


1.375(5) [27] 





Table 4 



Summary of infinite- volume estimates for the 3D antiferromagnetic RP 2 , 0(5) and 
0(4) models. We call rf to rf for RP 2 and to the rjr for O(N) models. FT stands 
for Field-Theory. 
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